Do current-density nonlinearities cut off the glass transition? 



o 
o 

> 

o 

:z; 

o 



I 

o 
o 



> 
o 



o 



X3 
O 

o 



X 



M. E. Cates^ and Sriram Ramaswamy 
^SUPA, School of Physics, University of Edinburgh, 
JCMB Kings Buildings, Edinburgh EH9 3JZ, United Kingdom 
^ Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, 
Bangalore 560012, India, and CMTU, JNCASR, Bangalore 560064, India 
(Dated: February 2, 2008) 

Extended mode coupling theories for dense fluids predict that nonlinear current-density couplings 
cut off the singular 'ideal glass transition', present in the standard mode coupling theory where 
such couplings are ignored. We suggest here that, rather than allowing for activated processes as 
sometimes supposed, contributions from current- density couplings are always negligible close to a 
glass transition. We discuss in schematic terms how activated processes can nonetheless cut off the 
transition, by causing the memory function to become linear in correlators at late times. 
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Mode Coupling Theory (MCT) offers an approximate 
but appealing account of the glass transition in simple liq- 
uids, whose physical content has long been debated. Ju- 
diciously used, particularly in colloidal liquids [1], MCT 
yields experimentally validated, semiquantitative predic- 
tions for the slowing down and then arrest of density fluc- 
tuations; yet when applied to much simpler systems, such 
as a single particle in an anharmonic trap, it yields wrong 
results [2] . It is thus important to delineate what the ap- 
proximations of MCT mean physically, so that these may 
be improved upon, and/or eventually unified with other 
theoretical approaches to glass formation. 

In this Letter we take this forward by examining the 
full or 'extended' versions of MCT (eMCT) in which non- 
linear couplings of the density field to kinetic degrees 
of freedom (currents) are considered. In two important 
and technically substantial papers. Das and Mazenko [3] 
and subsequently Goetze and Sjoegren [4] showed that, 
within both a field-theoretic and a projection-based for- 
mulation of eMCT (denoted eMCTl and eMCT2 respec- 
tively), such couplings cut off the dynamical singularity 
corresponding to the glass transition. (Some subsequent 
variants [5, 6], denoted eMCT3-5 are addressed below; 
an up-to-date review is provided by Das [7].) This cut- 
off is not present in the "standard" MCT (sMCT) which 
neglects nonlinear coupling to currents [8-11]. 

This cutoff was argued by some to represent an acti- 
vated hopping mechanism [12, 13]. The existence of such 
a mechanism is physically plausible: the activation en- 
ergy for relaxation should not exceed that required to 
expand a volume of order S,^ to a density low enough to 
melt the glass. Here ^ is the smallest length large enough 
that interactions across the boundary of such a melted 
droplet do not force it to refreeze in the same local state 
[14, 15]. Since sMCT predicts finite elastic moduli, this 
energy is finite, and the ideal glass transition is cut off, 
unless ^ is infinite throughout the glass phase. 

In recent years, however, the hopping interpretation 
of eMCT has seemed increasingly doubtful (see [7]). 



From a field-theory viewpoint MCT is formally a re- 
summed (one-loop self consistent) low temperature ex- 
pansion [3, 5, 16, 17]. In such an expansion, activated 
terms, which are exponential in /3 = l/ksT, represent es- 
sential singularities ( "instantons" ) and should not be ac- 
cessible even after resummation [18]. Second, the physics 
of eMCT should reduce to sMCT for interacting Brown- 
ian particles, which conserve neither momentum nor en- 
ergy [17, 19]. It would be surprising if there were no 
hopping processes in that case; yet MCT predicts none. 

Below we show that doubts about the physical basis 
of eMCT are well-founded. Accepting for these purposes 
that Refs.[3, 4] correctly implement their eMCTs, we ar- 
gue that the resulting cutoffs likely stem from new ap- 
proximations made in eMCT, beyond those of sMCT. We 
show that eMCT forms for the memory function 7fc(i), 
describing retarded friction, violate physical arguments 
based on a formal separation of current-density terms. 
We then discuss schematically how instantons [15, 20] 
may nonetheless cut off the sMCT glass transition. Our 
work may thus bring closer a unified description of col- 
lective arrest [13] and activated dynamics [14, 15, 21]. 

Our approach closely follows Zaccarelli et al. [22, 23]. 
We start with j — 1...N Newtonian particles of mass 
rn interacting via a nonsingular pair potential v(r) in 
volume V, at mean density po = N/V. The density is 
Pkit) = Ej e^''-'"^'^*^ and (1 - PoCk)-^ - Sk = {p-kPk)/N 
is the static structure factor. By standard arguments 
[22], an equilibrium memory function 7fe(t) exists so that 

Pk(t) + nlp^it) + Y^it) = fk{t) (1) 

where Ik = j^lk{t — t')p]i{t')dt'; /k is a noise of mean 
zero; and flj. = k"^ / (3mSk- Moreover, one may show from 
Newton's laws [22] that the memory function obeys 

7fe(i) = (/3m/7Vfc2)(^_k(0)(^kW + Y^{t))) (2) 

with J-k = ri^Pk + Tk+Qk; Tk = -E,(k.rj(t))2e"' '-^(*); 
and mFQk = - X^k' "^fe' C^^-^OPk-k' (i)Pk' (0- Suppose 
now that the particle velocities in T relax, on some finite 
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timcscalc t that docs not diverge near any glass transi- 
tion (below denoted 'rapid' relaxation), towards a global 
Boltzmann distribution, for which {riahf}) = ^aff/PnT- 
[22]. If so, then for t ^ t, Ti^{t) can safely be replaced by 
— Plait) /Pm. With this replacement (but noting that 
their Eq.ll fails without it) one recovers in combination 
Eqs.ll and 12 of [22]: that is, 7fe(t) = 7|(t) where 

7^(t) = {f3m/Nk^)[ik^poCk/Pmf{p-ia{0)pu{t)) (3) 
-(A;2poCfe//3m){(p_k(0)Qk(i)) + (Q-k(0)/3k(t)) 
+(p_k(0)yk(i))} + (Q-k(O)lk(i)) + (Q-k(O)Qk(t))] 

By treating momentum fluctuations as fast and inde- 
pendent of density, 7^ (t) neglects all the collective fluid 
modes responsible for current-density couplings; but, for 
t ^ T, nothing else of consequence seems neglected from 
(3) as derived in Refs.[22, 23]. As shown there, 7|(i) fur- 
ther reduces to the well known sMCT memory function 
Ik^'^'^it); if one treats Pk{t) as a Gaussian random pro- 
cess, and assumes that ( = (3vk + Ck = 0. We discuss such 
steps below, but without making either of them write 



7fe(i)=7fcW + A7fe(i) 



(4) 



where all "current-density couplings" reside solely in the 
second term. By this we mean that, although 7|(t) still 
formally depends on p (via Y{t)), it describes a system 
in which currents and densities are explicitly decoupled, 
in the sense that instantaneous equipa,rtition holds for all 
averages of equal-time products of particle velocity pairs 
in (2). The coupling term, Ajkit), whose explicit form 
follows by subtracting (3) from (2), then contains various 
correlators involving A7k = %^ + k^ p]^/ f3m. In real space 
this is AT(r,t) = V^V/s Y^ilriarlp-dapl l3m]5{Y-Yi{t)). 

Remarkably, both the eMCTl of [3] and the cMCT2 
of [4] yield approximations for the memory function that 
depart from the additive structure shown in (4). In the 
relevant regime (at low frequencies w, with k around the 
peak in Sk, and in the glassy regime where 7i<;(t) decays 
very slowly), the results are [3, 4, 24] 



, eMCTl 
Ik 



(co) = (l + zA,/u;)7^^^^(^) 



7^^^^(u.) 



(5) 
(6) 



Such results seemingly require that at low frequency the 
kinetic correction Ajk {^) either becomes larger than the 
standard contribution 7^(0;) (in the case of (5)), or re- 
mains comparable to it (in the case of (6)). We now argue 
that such behavior cannot be a physical consequence of 
the current-density couplings as defined above, because 
Ajkit) must remain negligible near a glass transition, 
even when the velocity and temperature fields of the fiuid 
(the only significant physics omitted from 7^(i)), are ac- 
counted for. 

To sec this, we define p{r,t),v(r,t) and T{r,t) = 
T + 6T as hydrodynamic number density, velocity and 



temperature fields, formed by local averaging over re- 
gions 6V{r) of volume containing i/(r, t) > I par- 
ticles [25]. The collective velocity field is v(r) = 
J2je6V(r) ^j/'^i^) ^^'^ density field is p{r) = j^(r)/A'^. 
Setting ri{t) —v{ri,t) = Svi{t), we next insist that relax- 
ation of the peculiar velocities 5vi (t) is rapid: local equi- 
librium alone allows the replacement m{Svia{t)Svifj{t)) ~ 
(1 — l/v[Y,t))5ai3kBT(jci,t) [26]. At temporal and spa- 
tial scales relevant to collective modes, we can then sub- 
stitute in (2) T(r,t) - V„V/3[p(r, t)w„(r, i)u;3(r, t) + 
{p{r,t) — X^^)kBT{r,t)/m]. This in turn can be writ- 
ten J2^=i Tii^^t), where mTi = mS7 ^{pVaVp); m72 = 
keTV'^p; and w.% = keiSTV^p + (p - \-^)V'^5T + 
(V5T)Vp]. These describe respectively gradients of hy- 
drodynamic inertial stress; gradients of ideal gas kinetic 
pressure; and a non-isothermal correction to the latter. 

We finally argue that heat conduction is rapid {5T ~ 0) 
in the relevant fc-range. This isothermal approximation 
is standard, even within eMCT [3-7], although it fails 
at lower k, where %>, causes sound-speed renormalization 
[27]. Using it, we find that AT(r,/) is merely Ti(r,t) = 
VaV/3/9WaW/3, SO that A7fc(t) contains only correlators 
of Ti(r,t) with V^p{r,t)] with Q(r,t); and with itself. 
The first of these terms, (V^p(0, 0)VaV/3pi'ai'/3(r, i)) 
in real space, cross-correlates gradients of density and 
inertial stress. But regardless of V^p(0, 0), the iner- 
tial stress term pVaV0{r,t) = J2i,je5V{r)'f'iarjf3/'^{r)>^^ 
rapidly decays (on time scale n, say) to its equipar- 
tition value kBTSa/s/mX'^; this is p-indcpcndcnt [25], 
hence uniform, and killed off by VaV/3. Accordingly 
(VV(O,0)Ti(r,t)) ~ for t » n. The only alternative 
is if relaxation of pv^Va is not rapid, but remains contin- 
gent on that of p itself, causing n to become very large in 
the glassy regime. Remarkably, a mode-coupling ansatz 
(VV(O,0)Ti(r,i)) ^ {W^p{Q,Q)V^p{v,t)){v^{v,f)) says 
exactly this. But such a result is wrong: it implies 
that, e.g., in window-glass, where atomic positions re- 
main unrelaxcd for centuries, the atomic momenta (of 
which V is a local average) fail to reach equipartition, 
even on that time scale. The same reasoning gives 
(Ti(O,0)Q(r,t)) 2± 0. It also gives (ri(0, 0)ri(r, t)) ~ 0, 
with one caveat; this contains a product of squares of two 
velocities at unequal times. It might therefore show alge- 
braic, rather than exponential, decay [28]. But such long- 
time-tails should play no role in the physics of glasses. 
(Indeed they are explicitly neglected in eMCT2 [4].) 

The above arguments complete our case that Ajkit) in 

(4), containing formally all current-density couplings, can 
be safely neglected anywhere near a glass transition. The 
results (5,6), which seemingly require Ajkit) to become 
large, miist then result from one of two things. Either 
these theories make injudicious approximations to Ajkit) 
and should be rejected; or they deal with the averages in 
(3) in a different manner from sMCT, so that (5,6) stem 
not from current-density couplings (as defined above) but 



3 



from different handling of the "standard" piece of the 
memory function, 7^(t). Recall that to reduce J^{t) to 
71*^*^-^(4), two further approximations are made in [22]. 
The first is to set C = 0, which amounts to adopting 
dynamically the Ramakrishnan Yussouff (RY) free en- 
ergy functional [7, 29]. The second is that Pkjio{t) may 
be " treated as gaussian" (a practice we below call "frag- 
ging" ) in the specific sense that its three-point correlators 
are set to zero, and its four-point correlators factored into 
sums of products of two-point correlators. A change to 
either assumption might give a cutoff, in principle. 

We argue next that eMCTl stems from injudicious ap- 
proximation to Ajk (t) rather than from an altered treat- 
ment of 7^(4). This is because eMCTl trags not only 
Pk, but also Vk and the momentum density gk: indeed it 
aims to determine selfconsistently the covariancc matrix 
of these tragged variables by an optimized resummation 
of nonlinearities. But if one believes RY is adequate [7] , 
and Ajk{t) is negligible, then the tragging of the density 
is sufficient to show via [22] that jk{t) = lk^^'^{t)- If so, 
then since the density remains tragged in eMCTl (trag- 
ging additional variables does not untrag existing ones), 
this theory can only give results different from sMCT by 
making additional approximations that sMCT avoids. 

Let us now consider how (5) actually arises. Within the 
formalism of Ref. [3] the renormalized equation of motions 
are GrJ V = fx, where the fluctuating (and tragged) 
variables y comprise p, v and g. Gxy is the propagator 
calculated in Ref. [3]. The hatted variables are auxiliary 
fields used to treat the noise terms fx (with variances set 
by G].l, and cross-correlations by G].^), and also to treat 
a constraint described below. Omitting noise terms for 
clarity, eMCTl 's equations of motion are [3] 

-iwpk(w) + ik.gk(w) = (7) 
-ia;gk(w) - ikc^pk(w) + 77fe(w)fc^Vk(w) = (8) 
gk(w) - /9oVk(w) + iQ!kpk(w) = (9) 

where a = A]./k'^ is a constant, and we have set m = 1 
so that p is now the mass density. These are linearized 
effective equations arising from a complicated nonlinear 
model; nonetheless one would normally expect them to 
have the usual hydrodynamic form, with renormalized 
coefficients. In this sense, Eqs. (7,8) are the expected 
equations of motion (with c a sound speed, which is k- 
independent in eMCTl); but Eq.(9) is not as expected. 
To understand it, note that in Ref. [3], a field v is deployed 
to maintain, as a nonlinear constraint, the defining equa- 
tion g(r, t) = p{r, i)v(r, t) of the fiuid velocity v. Despite 
this, the bare constraint g = po'v, instead of renormaliz- 
ing into the nonlinear one, g = pv, becomes (9), which 
in real space reads g = poV — aVp. Combining with (7) 
(i.e. p = — V.g) we obtain 

p = -PoV.v aV^ (10) 

instead of the continuity equation p = — V.(pv). Elimi- 



nating g, V from (7-9), one obtains 

-a;^Pk + c^Pk-«(wA;^/Po)%(w)[l + i^fe/a;]pk = /k (H) 

which confirms (5). Note we find no algebraic; error in 
Ref. [3] . The fault appears intrinsic to the whole program 
of extending MCT, which neglects vertex corrections [17], 
to currents. Here, instead of giving the vertex correction 
required to convert g = poV into g = pv, eMCTl creates 
a new (self-energy) term in a which, at least at the level of 
effective equations of motion, violates the constraint that 
V was supposed to enforce. The result (9) thereby intro- 
duces a parallel diffusion channel which requires, physi- 
cally, a momentum sink. While it might be tempting to 
associate the resulting 'motion' with activated hopping 
[4] , the end result (5) of eMCTl is incompatible with the 
physical content of (4) above, and we reject it. 

An alternative to eMCTl (call this eMCT3) in Ref. [5] 
likewise violates continuity; the authors 'justify' this on 
the grounds that eMCT3 also violates Galilean invari- 
ance. (An adequate theory should do neither.) The cut- 
offs predicted by eMCTl and eMCT3 are almost identical 
[7, 30]. Ref. [6] offers two further variants. In eMCT4, the 
constraint differs by a Jacobian from eMCTl; the results 
are unaffected. In eMCT5, the constraint is replaced by a 
scries expansion in Sp{r). The resulting equation of mo- 
tion for g is that found by substituting (9) in (8). Thus, 
so long as (8) is true, (9) is implied by eMCT5: the ex- 
pansion in Sp breaks the constraint g = pv in exactly the 
same sense as eMCTl, with the same consequences. 

We now turn to eMCT2. Here tragged variables are 
less easily identified [4, 31, 32], but an authoritative re- 
view [33] seems to suggest that the phase space density 
Pk{p,t) (with p particle momentum) is tragged. If so, 
since Pk(4) = X]pPk(p,t), and any sum of tragged vari- 
ables is de facto tragged, pk is tragged also in eMCT2. If 
this is true, then the cutoff (6) must either: (a) like (5), 
stem from injudicious treatment of A7fc(f) that sMCT 
bypasses; or (b) stem from avoiding use of = within 
the standard part of the memory function, 7^(4). In case 
(b), eMCT2 should be equivalent to Eq.l3 of [22], found 
by tragging the density in (3): this looks highly unlikely. 
The only remaining case is (c): that tragging of the den- 
sity is avoided, or partially avoided, in eMCT2. Result 
(6) might then remain a physical one; but, at least by our 
definition, one unrelated to current-density couplings. 

Our own view is that Ajk(t) is always negligible; RY 
(hence C = 0) generally sufficient; and the cutoffs in all 
forms of eMCT likely approximation-induced. Suppose 
one accepts this case: by what physical mechanism can 
the glass transition instead be cut off? Our results show 
that any cutoff mechanism must specifically invalidate the 
tragging of the density. We argue now that instanton-like 
activated processes ('hopping') do exactly this. 

When hopping occurs in a system that would other- 
wise be fully arrested, the generic relaxation process seen 
by any local neighborhood is as follows. First, wait a 
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long, stochastic time interval; then a locally large and 
rapid change occurs which reconfigures the density field 
nonperturbatively (an instanton). This is very unlike a 
Gaussian process, in which relaxations occurs in tempo- 
rally uniform small increments; tragging is appropriate 
only for relaxation by such gradual events. In the limit 
of pure instanton relaxation, the product of two fliictu- 
ating densities, p(0, 0)p(r, t), becomes suddenly random- 
ized after a long period of being constant; and crucially, 
/9^(0, 0)/9^(r, i) dccorrclatcs at exactly the same instant. 
Under these conditions, the (normalized) four-point cor- 
relator (p^(0,0)/9^(r,f)), which is almost but not quite 
the one needed in (3), is equal to the two-point correla- 
tor, not its square, in maximal violation of tragging. 

This suggests a schematic route by which sMCT might 
be cut off by hopping. (For more complete approaches 
along similar lines, see [15, 18, 20].) For t <^ tj (an in- 
stanton time) the minimal sMCT schematic model, called 
Fi2, is applied. This has a single correlator <I>(t) with 
memory function m{t) = Ai$(t) -|- A2$(t)^ and Ai,2 > 0. 
On increasing Ai this has [4], for A2 < 1, an unrealistic 
'type A' transition at AJ = 1; for 1 < A2 < 4, a realistic 

1 /2 

transition at A^ = 2A2 — A2; and for A2 > 4, the system 
is always in the glass. Suppose we generalise this to 

m{t) = Ai(u)$(t) + A2(u)$(t)2 (12) 

with u = t/rj. This defines parametrically some curve 
on the Ai,A2 plane; if this starts in the glass phase and 
ends in the fiuid, one can expect behavior appropriate 
to a cut-off glass (with details dependent on the curve 
chosen, and on the form of Ai(u)). However, treating the 
argument of the preceding paragraph literally would lead 
not only to A2(oo) =0, but also to Ai(oo) = Ai(0)-|-A2(0); 
this cannot give a trajectory of the required character [4]. 
This objection is perhaps not fatal since, in the 'scheme' 
that links F12 with fully k-dependent sMCT, the A's are 
proxies for certain k-intcgrals, whose crossover algebra in 
the instanton regime need not be as simple as assumed 
above. We leave detailed study of (12), and its possible 
relation to non-schematic models, for future work. 

On nearing completion of this work, a preprint [34] 
appeared giving more grounds to doubt the eMCT cutofF. 
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